THE ROLE OF INCIDENT PARTON TRANSVERSE MOMENTA IN HEAVY 

QUARK HADROPRODUCTION 



Yu.M.Shabelski 

Petersburg Nuclear Physics Institute, 
Gatchina, St. Petersburg 188350 Russia 

Talk, given at HERA Monte Carlo Workshop, 
"q^"' 27-30 April 1998, DESY, Hamburg 

a^ : 

. Abstract 

\ The conventional NLO parton model is enough for the description of total cross sec- 

^ ^' tions and one-particle distributions. In the case of two-particle correlations, the collinear 

■ approximation has failed, and it is necessary to account for the transverse momenta of 

^ \ initial partons. The different possibilities to do this are discussed. 
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^ ■ 1 Introduction 

^ : 

^ \ The investigation of heavy quarks production in high energy hadron collisions provides a method 

Q\ ' for studying the internal structure of hadrons. Some problems are the same in hadroproduction 

^ ■ and photo-/electroproduction processes. So, the review of the situation of the heavy quark 

hadroproduction can be useful for the interpretation of HERA data. 
Qh! In this talk we present a short review of heavy quark hadroproduction. The theoretical 

r-| ! predictions are usually obtained in the NLO parton model IQ. The assumptions which are used 

I for simplifications of the computations are considered in Sect. 2. In the case of one-particle 

distributions even LO (~ a^) parton model with collinear approximation is enough for the data 
description, NLO contributions (~ a^) only change the normalizations. On the other hand, in 
5^ ! the case of two-particle distributions, see Sect. 3, the collinear approximation has failed, and 

it is necessary to account for the transverse momenta of the incident partons. The possibility 
to include the transverse momenta of the incident partons in the framework of semihard theory 
[0, where the virtualities and polarizations of the gluons are taken into account, is considered in 
Sect. 4. In Ref. we presented the results for main and simplest subprocess, gg QQ (~ a^) 
for hadroproduction, and •jg —>■ QQ (~ a^) for photo- and electroproduction. 

2 Conventional NLO parton model 

The conventional NLO parton model expression for the heavy quark hadroproduction cross 
sections has the factorization form 

a{ab QQ) = Y1 j dxidxjGa/i{xi, fiF)Gb/j{xj, ij,F)a{ij ^ QQ) , (1) 



X 



where Ga/i{xi, hf) and Gb/j{xj, hf) are the structure functions of partons i and j in the colhding 
hadrons a and b, ^f is the factorization scale (i.e. virtualities of incident partons) and a{ij — 
QQ) is the cross section of the subprocess which is calculated in perturbative QCD. The last 
cross section can be written as a sum of LO and NLO contributions, d-{ij —>■ QQ) = (yU/j)cr^°'' + 
al{nii)alj\ where fiR is the renormalization scale, and cr|°^ as well as a^j^ depend practically 

only on one variable p = , s = XiXjSab- 

The expression (1) corresponds to the process shown schematically in Fig. 1 with 

qiT = g2T = . (2) 

The main contribution to the cross section at small x is known to come from gluon-gluon fusion, 
i= j = 9- 




Fig. 1. Heavy quark production in parton model. 

The principal uncertainties of any numerical QCD calculation of heavy flavour production 
are connected with the unknown values of the parameters: both scales, ^f and /i^Q, and the 
exact value of heavy quark mass, uiq. The values of both scales should be of the order of 
hardness of the considered process, however nobody can say what is better to use for scales, mg, 
rriT = ^ rriQ + or s. The phenomenological parton densities are sometimes (at very small 
x) in contradiction with the general properties of perturbative QCD. However it is just the 
region that dominates in the heavy quark production at high energies^. Another problem of 
parton model is the collinear approximation. The transverse momenta of the incident partons, 

^ These uncertainties should disappear when one sums up aU the high order contributions. Sometimes people 
say that strong scale dependence of the calculated results in LO or NLO means the large contribution of high 
order diagrams and weak dependence means their small contribution. Of course, it is not true. Strong scale 
dependence of NLO results means only strong scale dependence of high order contributions but at some fixed 
scale value the last ones can be numerically small. Weak scale dependence of NLO results means weak scale 
dependence of high order terms but they can be numerically large. 

^In the case of charm production, = 1.4GeV, at LHC, = 14 TeV, the product 2:1X2 of two gluons (both 
Xi and X2 are the integral variable) is equal to 4 • 10~*. 
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qiT and qjT are assumed to be zero, and their virtualities are accounted for only via structure 
functions; the cross sections a\°^ and alj'' are assumed to be independent on these virtuahties. 

The NLO parton model calculations of the total cross sections of cc and bb production, as 
functions of the beam energy, for Tt~N and p — N collisions can be found in 0. These results 
depend strongly (on the level of several times) on the numerical values of quark masses as well 
as on the both scales, fip and hr. Some experimental data are in contradiction with each other, 
however generally they are in agreement with NLO parton model predictions. 
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Fig. 2. pT^-distributions for pp ^ b + X at y/s = 1.8 TeV, for different values of rapidity. 



The NLO contributions to one-particle distributions lead only to renormalization of LO re- 
sults, practically without correction of the shapes of a distributions |]. It means that instead 
of more complicate calculation of pt, or rapidity distributions, in NLO, it is enough to calculate 
them in LO, and multiply after by K-factor 

A-^^^^. (3) 

which can be taken, say, from the results for total production cross sections. The comparison of 
LO + NLO calculations with LO multiplied by K-factor is presented in Fig. 2 taken from Ref. 

0- 

The values of K-factors and their energy and scale dependences for several sets of structure 
functions were calculated in Refs. pUf . 



The experimental data for a;i?-distributions of D- mesons produced in nN interactions [|TT], |T2 
are in agreement with the parton model distributions for bare quarks, as one can see in Fig. 3 
taken from H]. It means that the fragmentation processes are not important here, or they are 
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Fig. 3. Experimental xp distributions for D mesons, compared to the NLO parton model 

prediction for charm quarks. 



compensated by, say, recombination processes. The shape of XiT'-distributions does not depend 
practically on the mass of c-quark. 

The data on one-particle pT-distributions, including the hadronic colliders data for the case 
of beauty production, also can be described by the NLO parton model, see 0. 

3 Azimuthal correlations and failure of the collinear ap- 
proximation 

The azimuthal angle is defined as an opening angle between two produced heavy quarks, 
projected onto a plane perpendicular to the beam. In the LO parton model this angle between 
them is exactly 180°. In the case of NLO parton model a distribution over angle appears 0. 

The investigation of such distributions is very important. In one-partricle distributions, the 
sum of LO and NLO contributions of the parton model practically coinsides with the LO con- 
tribution multiplied by i^-factor. So we can not control the magnitudes of LO and NLO contri- 
butions separately. In the case of azimuthal correlations all difference from the trivial — tt) 
distribution comes from NLO contribution. 



The experimental data on azimuthal correlations are claimed (see |0 and Refs. therein) 
to be in disagreement with the NLO predictions, for the cases of charm pair hadro- and pho- 
toproduction at fixed target energies. The level of disagreements can be seen in Fig. 4 (solid 
histograms) taken from Ref. 0. These data can be described |0, assuming the comparatively 
large intrinsic transverse momenta of incoming partons {kx kick). For each event, in the longi- 
tudinal centre-of-mass frame of the heavy quark pair, the QQ system is boosted to rest. Then 
a second transverse boost is performed, which gives the pair a transverse momentum equal to 
Pt{QQ) + ^t(I) + kTi^); kril) and /ct(2) are the transverse momenta of the incoming partons, 
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Solid : QCD NLO 

Dashed: NLO+k, kick, <k?> = 0,5 GeV^ 
Dotted: NLO+k, kick, <kT>=l GeV^ 
X : TA75 data (ttN) 
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Solid : QCD NLO 

Dashed: NLO+kT kick, <k?>=0.5 GeV'^ 
Dotted: NLO+kT kick, <k?> = l GeV'^ 
X : WA92 data (ttN) 
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Fig. 4. Azimuthal correlation for charm production in vriV collisions: NLO parton model and 
kx kick calculations versus the WA75 and WA92 data. 



which are chosen randomly, with their moduli distributed according to 

1 dN 



eM-KKK))- (4) 



The dashed and dotted histograms in Fig. 4 correspond to the NLO parton model prediction, 
supplemented with the kx kick, with {k"^) = 0.5 GeV^ and (/cf.) = 1 GeV^, respectively. We see 
that with (kx) = 1 GeV^ it is possible to describe the data. 

However the large intrinsic transverse momentum significantly changes one-particle p^-distri- 
butions of heavy fiavour hadrons, which were in good agreement with the data. The solid 
curves in Fig. 5 taken from [^] represent the NLO parton model predictions for charm quarks 
pr-distributions which are in agreement with the data. The effect of the kx kick results in a 
hardening of the spectrum. On the other hand, by combining the fcy kick with (k^) = 1 
GeV^ and the Peterson fragmentation [Q, the theoretical predictions slightly undershoot the 



data (dot-dashed curves) 

The kx kick can only very weakly change the Xi?-distributions of produced c-quarks. Fig. 3, 
and after accounting the fragmentation these distributions can become too soft. 

Let us consider why the conventional NLO parton model with collinear approximation works 
reasonably for one-particle diustributions, and, at the same time, it is in disagreement with the 
data on azimuthal correlations. 

The contribution of the processes of Fig. 1, which governs the heavy quark production can 
be written^ as a convolution of initial transverse momenta distributions, I{qiT) and I{q2T), with 
squared modulo of perturbative QCD matrix element, |M(gir, g2T,PiT,P2r)P : 

c^qcd{QQ) oc j d^qiTd^q2Tl{qiT)I{<l2T)\M{qiT,q2T,PiT,P2T)? ■ (5) 



^We omit for simplicity all factors which are non-essential here. 
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Fig. 5. Charm distribution measured by WA92 and E769, compared to the NLO parton 
model predictions, with and without the non-perturbative effects. 



Now there are two possibihties: 

i) the essential values of initial transverse momenta are much smaller than the transverse mo- 
menta of produced heavy quarks, qiT -C Pix, and 

ii) all transverse momenta are of the same order, ~ piT- 

In LO for the first case (i) we can wait that one-dimentional pt distributions should be 
more broad than the distributions on the transverse momenta of the quark pair, because in LO 
QiT + q2T = Pit + P2T- NLO gives here a correction, numerically not very large, see Fig. 2. In 
this case one can replace both the initial distributions /(^it) by 5-functions, This reduces 

the expression (5) to the very simplified one : 

(TcollXQQ) « |M(0,0,PiT,P2T)r , (6) 

in total agreement with Weizsaecker- Williams approximation in QED. 

In the second case (ii) the distributions on the transverse momenta of the quark pair should be 
the same, or even more broad than the one-particle p^-distributions, and namely this situation is 
realised ^ . In this case we can not wait a priory that the Weizsaecker- Williams approximation 
will give good results, however it works quite reasonably in the case of one-particle distributions. 
In the case of distribution on the transverse momentum of the heavy-quark pair, p%{QQ) we 
measure (only approximately in NLO) the distribution over the sum of transverse momenta of 
incident gluons. In this case it seems to be senseless to replace the distributions /(q'it) and 
I{(12t) by 5-functions, and to expect a reasonable agreement with the data. The same can be 
said about the azimuthal correlations. 

The kt kick |^ effectively accounts for the transverse momenta of incident partons. It uses 
the expression which can be written symboUically as 

(THckiQQ) oc I{qiT)I{q2T) ® |M(0,0,piT,P2T)r , (7) 

and the main difference from the general QCD expression Eq. (5) is that due to absence of q^T 
in the matrix element the values of (fc^) in Eq. (4) should be different for different processes and 
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kinematical regions. The reason is that in Eq. (5) the values I{qiT) decrease at large Q^jfj^ as a 
weak power (see next Sect.), i.e. comparatively slowly, and more important is the g|p dependence 
of the matrix element. In the last one the corrections of the order of Q'^rZ/i^, where //^ is the 
QCD scale, are small enough when ^'j^'//^^ << 1 and they start to suppress a matrix element 
value when q^rp/ ~ 1. 



4 Heavy quark production in semihard approximation 

Let us consider another approach, when the transverse momenta of incident gluons in the small-x 
region appear from the diffusion of transverse momenta in the gluon evolution^. This diffusion 
is described by the function ip{x, q^) determined p as 



V9(x, g^) = 4A/2 7r 



d [xG{x, q^ 



dq^ 

where G(x, q^) is usual gluon structure function. 

In principle the function ip{x, q^) which determine the probability to find gluon with fixed 
value of transverse momentum, g^^, depends on three variables, x, qx and gluon virtuality q^. 
However at small x in LLA q^ ^ ~q^i and it leads to comparatively weak dependence of (^{x, q^) 
on q^ (strongly different from exponential dependence in Eq. (4)) due to weak g^-dependences 
of phenomenological structure functions. 

The exact expression for gluon gj^-distributions can be obtained, as a solution of the nonlinear 
evolution equation. The calculations |jl8[ result in difference from our ip{x, g^) function only 
about 10-15%. 

The matrix element Mqq accounting for the gluon virtualities and polarizations is much more 
complicate than the parton model one. That is why we consider only LO contribution of the 
subprocess gg — >■ QQ. The differential cross section of heavy quarks hadropro duct ion has the 
form 

j d^qiTd^q2T5{qiT + q2T-ViT-P2T) 



dcTpp _ 1 1 / 72 ^2 



dytdy^d'^PiTd'^P2T i^n)^ (s)^ 

X — 2 —v{<li,ym<l2,x)WQQ\ . (9) 

"?! Q'2 

Here s = IpaPb and yl 2 are the quarks' rapidities in the hadron-hadron c.m.s. frame. 

Eq. (9) enables to calculate straightforwardly all distributions concerning heavy flavour one- 
particle, or pair production. However there exists a problem coming from infrared region. Gluon 
structure function in Eq. (8) is not determined at small virtualities, so the function ip{x, g|) is 
unknown at the small values of g| and q\. To solve this problem we will use the direct consequence 
of Eq. (8) [0 

1 n"" 

xG{x, g2) = xG{x, Ql) + -j=— / dqjip{x, ql) , (10) 



^0 



*The similar approach based on fcr-factorization formulae can be found, for example, in Refs. Eq, Eq, p 
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and rewrite |Q the integrals in the Eq. (9) as the sum of four contributions. The first one is 
determined by the product of two gluon distributions, G{x, Ql) and Giy, Qo), and it is the same 
as the conventional LO parton model expression. Next three terms contain the corrections to 
the parton model. If the initial energy is not high enough, the first term dominates. In the case 
of very high energy the first term can be considered as a small corrections, and our results are 
differ from the conventional ones. In the cases when the collinear approximation is available, 
our results only slightly differ from the parton model. 




5 10 15 20 25 30 35 40 



Fig. 6. Cross section of beauty production in CDF. 



It is illustrated in Fig. 6 taken from 0, where we compare our calculations of 6-quark p^- 
distributions with the experimental results of CDF collaboration and with two variants of parton 
model calculations. The distributions over the azimuthal angle can be found in ||2^ . 



The essential values of qiT and g2T in our calculations increase with increase the value of p^*" 
of detected 6-quark. In the language of kick it means that the values of {k'^) will be also 
increased. 



5 Conclusion 

The experimental results on total cross sections for charm and beauty production are in agree- 
ment with the conventional parton model predictions, using reasonable values of QCD scales 
and quark masses. The data ow xp and pt distributions are also in reasonable agreement with 
parton model without any fragmentation functions^. 

^The Pt distributions of charm photoproduction measured by E691 Coll. are more hard than the NLO 
parton model predictions that can be considered as an argument for including a non-perturbative 
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Moreover, the shapes of one-particle LO and NLO distributions practically coinside. It means 
that instead of calculation the NLO contributions, it is enough to calculate only LO contribu- 
tions, and rescale them using K-factor taken, say, from the calculated ratio of total cross sections. 

In the case of distribution over the total transverse momentum of the produced quark pair, 
or azimuthal correlations, the conventional NLO parton model with collinear approximation can 
not describe the data. The kx kick ^ allows one to describe these data, however the problems 
with one-particle pT^-distributions appear, which can be solved by introducing the fragmentation 
function. The last way should produce the problems in description of a;ir-spectra. Moreover, it 
seems that the (k^) values should depend on the process and the kinematical regions. 

Another possibility to solve the problems of initial transverse momenta is to use semihard 
theory, accounting for the virtual nature of the interacting gluons, as well as their transverse 
motion and different polarizations. It results in a qualitative differences with the LO parton 
model predictions In Ref. the values of F2{x, Q"^) were calculated using phenomeno- 



logical gluon structure functions, and the infrared contributions to F2{x,Q'^) were investigated 
in details. The possible estimations of the shadow corrections in the processes of heavy flavour 
production can be found in Refs. p, p4[ . 

I am grateful to M.G.Ryskin and A.G.Shuvaev for multiple discussions, to M.L.Mangano for 
very useful critical comments and to E.M. Levin who participated at the early stage of this ac- 
tivity. I thank the Organizing Commitee of HERA Monte Carlo Workshop for financial support. 
This work is supported by grant NATO OUTR.LG 971390. 
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